An integral equation approach to effective interactions between polymers in solution 
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We use the thread model for linear chains of interacting monomers, and the "polymer reference 
interaction site model" (PRISM) formalism to determine the monomer-monomer pair correlation 
function h mm (r) for dilute and semi-dilute polymer solutions, over a range of temperatures from 
very high (where the chains behave as self-avoiding walks) to below the 9 temperature, where phase 
separation sets in. An inversion procedure, based on the HNC integral equation, is used to extract 
the effective pair potential between "average" monomers on different chains. An accurate relation 
between /i mm (r), h cc (r) [the pair correlation function between the polymer centers of mass (cm.)], 
and the intramolecular form factors is then used to determine h cc (r), and subsequently extract 
the effective c.m.-c.m. pair potential v cc (r) by a similar inversion procedure. v cc (r) depends on 
temperature and polymer concentration, and the predicted variations are in reasonable agreement 
with recent simulation data, except at very high temperatures, and below the 9 temperature. 



I. INTRODUCTION 

The Statistical Mechanics of dilute and semi-dilute so- 
lutions of interacting polymers may be greatly simpli- 
fied by adopting a coarse-graining strategy for this com- 
plex many-body system, whereby the initial monomer or 
segment-level description is reduced to a "soft-colloid" 
representation involving only the centers of mass (cm.) 
of the polymer coils. This is formally achieved by trac- 
ing out the individual monomer degrees of freedom for 
fixed positions of the c.m.'s of the polymers. This reduc- 
tion procedure results in effective interactions between 
the c.m.'s which are state dependent, since they are in 
fact free energies of the "fluid" of monomers, for each 
configuration of the c.m.'s of the N interacting coils. In 
the low concentration limit, this free energy reduces to a 
pairwisc additive sum of effective potentials of mean force 
acting between the c.m.'s of overlapping polymers, which 
depend only on temperature. Such pairpotentials were 
first considered by Flory and Krigbaum p| , and later put 
on a firmer basis by Grosberg et al. Q , who showed that 
the effective pair potential between the c.m.'s of linear 
self-avoiding walk (SAW) polymers remains finite at full 
overlap in the scaling limit (where the number of seg- 
ments L — > oo), and of the order of a few /c^T, reflecting 
the purely entropic origin of the effective interaction in 
the athermal SAW model. Quantitative estimates of the 
effective cm. pair potential between two isolated SAW 
polymers were obtained by Monte Carlo CMC) simula- 
tions of on- and off-lattice models S0, M E3, and 
by renormalization group calculations 8]. These studies 
show that the zero concentration pair potential is purely 
repulsive and reasonably well represented by a Gaussian 
of amplitude ~ 2fcsT and width equal to the radius of 
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gyration R g . 

The SAW model describes polymers in good solvent. 
Poor solvent conditions may be modeled by adding a fi- 
nite attractive interaction between monomers and lower- 
ing the temperature. MC simulations jj,19| show that the 
resulting effective cm. pair potential becomes less repul- 
sive as the temperature T is lowered, and can become 
attractive in the vicinity of the temperature. 

When the polymer concentration is increased, three- 
and more-body effective interactions between the poly- 
mer c.m.'s come into play. These can also in principle 
be calculated by MC simulations of clusters of three or 
more polymers |Io|. but this systematic procedure be- 
comes rapidly untractable, since the effective interactions 
depend on an increasing number of variables. 

An alternative approach is to restrict the effective in- 
teractions between polymer c.m.'s at finite concentra- 
tions to a sum of pair potentials, the latter being deter- 
mined by an inversion procedure of the cm. pair distri- 
bution function based on integral equations for the pair 
structure of classical fluids 0> B E} • The price to pay 
is that the resulting pair potentials now depend on poly- 
mer concentration, in addition to temperature. The vari- 
ation with concentration is moderate, but significant, in 
the case of the SAW model for polymers in good solvent 
0- 0] ,and becomes very strong for polymers in poor 
solvent 

Since the effective pair potentials are obtained by in- 
verting cm. pair correlation data from MC simulations, 
this strong variation with polymer concentration means 
that full monomer-level simulations must be carried out 
for each thermodynamic state of the polymer solution, 
thus partly defeating the purpose of the coarse-graining 
strategy. 

It is hence important to explore numerically less de- 
manding theoretical methods for the determination of 
the required polymer pair correlation functions. One ob- 
vious candidate is the polymer reference interaction site 
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model (PRISM) theory, which has proved very success- 
ful to determine the pair structure of interacting poly- 
mer solutions and melts E2- PRISM grew out of the 
RISM integral equation of Chandler and Andersen 12] 
for the pair structure of simple molecular fluids. The the- 
ory includes intramolecular connectivity via the macro- 
molecular form factor, and provides accurate monomer- 
monomer pair distribution functions, g mm (?). The c.m.- 
c.m. pair distribution function, g cc (r), which needs to be 
inverted to extract the effective cm. -cm. pair potential, 
v C c( r ), can be derived from g mm (r) by a simple, accu- 
rate relation which was established recently within the 
PRISM framework 13]. 

This strategy, which is highly economical on compu- 
tational resources, is explicitly carried out in this paper, 
on the basis of the Gaussian thread model for interact- 
ing polymers El EH EE El- PRISM theory, in con- 
junction with an approximate closure relation, allows the 
monomer- monomer pair distribution function g mm (r) to 
be calculated analytically or numerically depending on 
the expression for the polymer form factor. The cm. pair 
distribution function is then extracted from g m m{f) using 
the relation of Ref. [13j together with appropriate form 
factors, and v cc (r) is finally obtained by hypernetted- 
chain (HNC) inversion [a LA El- The variation of the 
effective pair potential with temperature and polymer 
concentration is then examined in detail, and compared 
to the "exact" results extracted from MC simulations 

siaa. 



II. POLYMER MODEL AND PRISM THEORY 

We consider solutions of N linear polymers made up of 
L freely jointed monomers or segments of diameter d and 
spacing a. Monomers belonging to the same or different 
chains are not allowed to overlap, i.e. they experience 
an infinite mutual repulsion for monomer-monomer dis- 
tances r < d, and, for r > d, they attract each other via 
a Yukawa interaction, 



v(r) 



-e—e 
r 



(1) 



where e > is the energy scale and a the range of 
the monomer-monomer attraction [TH Il6| . The rele- 
vant thermodynamic variables are the temperature T [as 
usual, we define (3 — (fc^T) -1 ] and the polymer concen- 
tration p (number of polymer coils per unit volume); the 
corresponding monomer density is p m = pL. As usual, 
the spatial extension of the chain will be characterized by 
its radius of gyration R g and the polymer overlap con- 
centration, at which different chains start to overlap, is 
defined as p* — (^irRg)^ 1 . In this model of polymer so- 
lutions, the attractive interaction (|TJ represents solvent 
mediated monomer-monomer interactions in an effective 
way. Then, with respect to intermolccular interactions, 
large T values correspond to good solvent conditions, 
while lower T values mimic poor solvent conditions. 



PRISM theory assumes all monomers to be equiva- 
lent regardless of their position along the polymer chain. 
Its aim is to calculate the site-averaged pair distribu- 
tion function g mm (r) of monomers on different chains. 
To this end one introduces the pair correlation func- 
tion, h mm {r) = g m m{f) — 1, the direct correlation func- 
tion c mm (r), and their Fourier transforms, respectively 
hmm{q) and c mm (q), which are related by the follow- 
ing generalization of the Ornstcin-Zernike (OZ) equation 

llj rnrn 

(?) 

Program 

(?)]. (2) 

where u) mm (q) is the form factor of the polymer chains. In 
the following, it will be convenient to introduce the aux- 
iliary quantities tl mm (q) = Q mm {q)/L [hence, Sl mm {q = 
0) = 1] and C mm (q) = L 2 c mm (q), such that Eq. J2J may 
be rewritten as 

h mm (q) = ^ mm (q)C mm (q)[Ci m m(q) + ph mm {q)\- (3) 

The OZ equation, Eq. must be supplemented by a 
closure relation for C mm (q), as well as by a prescription 
for VL mm (q). While self-consistent schemes exist to com- 
pute the latter as a function of temperature and density, 
they are quite complex to implement |ll| . In the present 
work, we will make the simplifying assumption that the 
polymers obey Gaussian statistics. This approximation, 
which is in principle only valid in the melt or close to the 
9 temperature, is mainly needed for technical reasons, 
due to the lack of analytical results for non-Gaussian 
polymers. Thus, the form factor O mm (q) is given by 
the well-known Debye ex pre ssion for ideal chains in the 
continuum limit, namely |18| 
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Replacing this exact form by its simple Pade approxi- 
mant, 



^mm (?) 



1 



1 + x 2 /2 



, x = qR. 



g> 



(5) 



allows analytic results to be obtained within the thread 
limit [lllialig, but we will see in the following that it 
can lead to significant problems in the present context. 

In this work, we are interested in a mesoscopic level of 
description of polymer systems. For simplicity, it is thus 
useful to get rid of the microscopic details of the polymer 
model. A way of doing this in the PRISM framework is by 
taking the so-called thread limit. In the case of excluded 
volume interactions only, it is defined by letting d and a 
go to zero and L to infinity, keeping the ratio T = a/d , 



ig the ratio 1 = i 

El 13 IB E 



R g and p/p * = ^irpR? g constant 

In Refs. E3 E| j the attractive part of the monomer- 
monomer interaction potential was left untouched when 
taking the thread limit. However, for simple polymers, 
it turns out that the range of this attractive tail scales 
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FIG. 1: Density dependence of the direct correlation function 
Co, in the presence of excluded volume interactions only. For 
the exact form factor Q, Co has to be numerically determined, 
while for the approximation JSJ, the analytic result IA2I is 
plotted. 



with the monomer size. Thus for this work it seemed 
physically appropriate to supplement the previous lim- 
iting process with the limit of a vanishing range of the 
attractive potential, namely a — > 0, keeping the ratio 
A = a/a constant. This can be interpreted as a special 
limit of the theory of Refs. EEl, as shown in Appendix 

To escape the trivial case of non-interacting polymers, 
the thread limit requires a proper formulation of the clo- 
sures of the PRISM-OZ equation as well. In this work, 
we have applied the so-called RMPY/HTA closure, which 
has been found to be the most satisfactory in earlier work 
UHdUm The details are given in Appendix lAl and we 
only quote here the equations relevant to our calculation. 
The monomer-monomer correlations are thus given by 



h>mm (q) — 



cR 3 g Q mm (q) 2 
1 - pcR 3 Sl mm (q) 



with 



C = Co 1 



(6) 



(7) 



where is the 9 temperature of the polymer model, cq, 
which depends on pR? g only, is obtained by imposing that, 
in the presence of excluded volume interactions only (i.e. 
for infinite temperature), h mm obeys the hard core con- 
dition, h mm (r = 0) = — 1. It is plotted for both chain 
form factors, Eq. 10} and (J5J, in Fig. [J 

With the above monomer pair correlation functions, ef- 
fective pair potentials between polymer chains may now 
be determined. We first consider in the next Section 



the situation where the interaction site for the effec- 
tive interaction between polymer coils is chosen to be 
a monomer; subsequently effective interactions between 
polymer c.m.'s are discussed in Sec. II VI 



III. EFFECTIVE INTERACTION BETWEEN 
MONOMERS 

Starting from the original polymer system, the choice 
of the effective coordinates used to specify the position 
of each polymer coil, and relative to which the coarse- 
graining procedure is performed, is largely arbitrary. 
Most of the time, the center of mass is considered, but it 
has sometimes been suggested to use specific monomers, 
like the midpoint or the endpoints, as in Refs. [SB 112 • 
This is a very natural choice for star polymers as well, 
where the star center provides an obvious reference point 
[23|. In this Section, we thus attempt to build such a 
monomer-monomer (or site-site) effective interaction. 

In the present PRISM-based approach, by construc- 
tion, no distinction is made between the different 
monomers on a chain. However, one can build an ef- 
fective monomer-monomer interaction based on the con- 
cept of an average monomer [24| • The idea goes as fol- 
lows. On each chain, a random monomer is tagged, and a 
coarse-grained polymer-polymer interaction is computed 
by tracing out, for fixed (spatial) positions of the tagged 
monomers, all the other monomer degrees of freedom. 
Averaging this free energy over the random locations of 
the tagged monomers on their chains then leads to a 
quantity which can be given the meaning of an effective 
interaction between average monomers. 

In practice, an effective pair potential v aa (r) may be 
extracted by an inversion procedure, based on integral 
equations for the pair structure of simple fluids (2^| . For 
this purpose, we consider a fluid of average monomers, 
of density p (each polymer chain is originally associated 
with a single tagged interaction site), the pair correla- 
tions of which are characterized by h aa (r) = h mm (r), 
where h mm (r) is the monomer-averaged correlation func- 
tion determined in the previous section [2g. A direct 
correlation function c aa (r) may be extracted from h aa (r) 
via the OZ relation for simple (atomic) fluids which in 
Fourier space reads 



h aa (q) = c aa (g)[l + ph aa (q)]. 



(8) 



The effective interaction potential may finally be deter- 
mined using one of the standard closures for simple flu- 
ids. Since the effective interactions are expected to be 
relatively weak (in particular, they contain no hard core 
repulsion), we follow earlier work IS II II 13, and adopt 
the very accurate HNC closure, which leads to the fol- 
lowing expression for v aa (r): 

f3v aa {r) = h aa (r) - c aa {r) - ln[l + h aa (r)]. (9) 

With this closure, a fully analytic expression for f3v aa (r) 
can be obtained if the approximate expression JSJ for the 
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chain form factor is used. Details are given in Appendix 

m 

We first consider the case of excluded volume interac- 
tions only, i.e. c= cq. The analytic Pv aa (r) is plotted in 
Fig. [3] for various polymer densities. It diverges at zero 
separation and, when the density increases, this diver- 
gence manifests itself at shorter distances and becomes 
steeper. At the same time, the repulsive tail of the poten- 
tial becomes longer-ranged, as can be clearly seen on the 
plot of 47r(r / R g ) 2 /3v aa (r), a quantity playing a crucial 
role for the calculation of the thermodynamic properties 
of fluid systems characterized by similar "soft" interac- 
tion potentials [27j • These results are in good qualitative 
agreement with the simulation data of Bolhuis and Louis 
for midpoint and endpoint effective pair potentials Q . 

Recently, much attention has been devoted to the de- 
velopment of analytic expressions for effective pair po- 
tentials between the centers of star polymers in solution 
|23|. Since a linear polymer, seen from one of its end- 
points or from its midpoint, can be identified with a one- 
or two-arm star, respectively, it is interesting to compare 
those calculations with the present theory. 

Since, according to Eq. i|A3|l . 1 + h mm (r) vanishes lin- 
early for r — > 0, [3v aa {r) diverges like — \n(r/R g ). Such a 
weak, logarithmic divergence is present in the star poly- 
mer models, where it has been derived using scaling ar- 
guments, but the prefactor, which should be given by 
the two-arm result in the limit of very long chains where 
end effects are expected to be negligible, appears overes- 
timated by the present theory compared to scaling argu- 
ments (which predict a prefactor equal to 5/ 3 / 2 /18 for 
a star with / arms) and renormalization group results 
|28|. This quantitative disagreement is not surprising, 
since the present calculation relies on the use of Gaus- 
sian form factors in good solvent conditions. However, a 
similar calculation done with an approximate form factor 
reproducing the correct asymptotic behavior of f2 mm (<z) 
in good solvent still does not provide the correct value 
and this failure was in fact argued to be inherent to the 
thread PRISM approach [l9| . 

For large distances, (3v aa (r) asymptotically takes a 
Yukawa form [2{J, 

/3v aa (r) cx — — . (10) 
r/Rg 

Such a long range behavior was postulated for star poly- 
mers, on the basis of the blob picture of the star in- 
ner structure. However, because the blob picture breaks 
down for small /, it was advocated that this asymptotic 
behavior should not be valid for small functionalities, 
/ < 10, and should be replaced by a Gaussian tail, by 
analogy with the behavior of the effective pair potential 
between chain centers of mass ■ The present calcula- 
tion does not appear to support this choice. 

Finally, and maybe more importantly, the present cal- 
culation reveals that Pv aa {r) changes significantly with 
density. For instance, the range of the asymptotic 
Yukawa tail more than doubles between p/p* = and 



p/p* = 1. If such a behavior also applies to star poly- 
mers, it could have a major impact on the calculation of 
the thermodynamic properties of these systems, in par- 
ticular on their phase diagrams. 

In Fig. the numerically determined (3v aa {r) for the 
exact Debye form factor is compared to the analytic 
result obtained with the Pade approximant JSJ, at var- 
ious polymer densities. A thorough analysis shows that 
the differences between both are only quantitative. Both 
share the same — \n(r/R g ) divergence for small separa- 
tions, and both display an asymptotic Yukawa tail, which 
is of shorter range when is used. 

We now turn to the effect of the attractive monomer- 
monomer interaction with changing temperature. In 
Fig. the analytic (3v aa (r) calculated with the ap- 
proximate chain form factor is plotted at density 
pRg = 0.07 for various values of the temperature. Sim- 
ilar results are obtained at all densities, and with the 
exact Debye form factor as well. The effective potential 
is seen to become less repulsive when the temperature is 
decreased from oo to O, it vanishes identically at O and 
grows more and more attractive if the temperature is fur- 
ther decreased. This trend reverts when the "spinodal" 
line is approached, leading to non-monotonous, mostly 
repulsive effective potentials of increasing range. 

But a word of caution is needed here. Indeed, as shown 
in Appendix^ at all finite temperatures, the hard core 
condition h mm (r = 0) = —1 is no longer obeyed. Thus, 
the effective pair potential between average monomers 
does not diverge at full overlap as it should, if this con- 
dition was consistently enforced at all state points. This 
might not be so important, since, as it has already been 
mentioned, the product Ait{r / R g ) 2 (3v aa {r) is in general 
more relevant than the potential itself, and in this quan- 
tity the short range details of the potential are strongly 
reduced by the vanishing prefactor. Another source of 
concern is that it is not known to what extent the re- 
sults obtained below when approaching the "spinodal" 
line are physically meaningful, since the corresponding 
remnant of a binodal line, which would properly iden- 
tify the domain of existence of the low density phase, is 
not known in the present theory [3(| ■ For these reasons, 
the present results are only given for reference and would 
need to be validated, both by implementing alternative 
integral equation schemes or by investigating /3v aa (r) as 
a function of temperature by molecular simulation. 

To summarize, we have shown in this Section that in- 
tegral equation theories for macromolecular systems, like 
the present thread PRISM theory, could be useful tools 
to compute monomer-based effective interaction poten- 
tials between polymers. At least in the case of excluded 
volume interactions only, they lead to very reasonable 
predictions for the effective interactions between average 
monomers on polymer chains, compared to the presently 
available simulation and theoretical results. 

But, the concept of an average monomer, which is 
needed in the present theory because, at variance with 
star polymers, there is no obvious choice of a refer- 
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FIG. 2: Effective pair potentials between average monomers (left-hand frame), obtained with the approximate polymer form 
factor ©, at various reduced densities pR 3 g , in the presence of excluded volume interactions only; the right-hand frame shows 
the potentials multiplied by An{r / R g ) 2 . 




FIG. 3: Comparison of the effective pair potentials be- 
tween average monomers [multiplied by An(r / R g ) 2 ] obtained 
with the exact and approximate polymer form factors, Q 
and © respectively, at various densities, in the presence 
of excluded volume interactions only. From left to right: 
pR 3 g = 0, 0.28, 0.56, 1.12. 



FIG. 4: Effective pair potentials between average monomers, 
obtained with the approximate polymer form factor 0, at 
various temperatures, for the reduced density pR g = 0.07. 



ence monomer, is rather artificial and in practice con- 
ceals significant technical difficulties. Indeed, this con- 
cept has a simple meaning only in the limit of very long 
chains, where all monomers can be considered as equiva- 
lent. In Ref. 0, the effective interactions between central 
monomers and between end monomers have been found 
to differ considerably, showing that strong finite length 



effects should be expected on effective interactions be- 
tween average monomers. Thus, for a better control of 
these finite size effects, an alternative approach seems de- 
sirable, and it has been found that considering effective 
interactions between polymer c.m.'s is appropriate in this 
respect 0, S H E . 
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IV. EFFECTIVE INTERACTION BETWEEN 
CENTERS OF MASS 

As has been shown in recent Monte Carlo studies [E 
IE Q an effective pair potential v cc (r) between the 
c.m.'s of polymer chains may be extracted by an inversion 
procedure analogous to the one described in the previous 
Section. To that purpose one considers a simple fluid of 
density p, the pair correlations of which are characterized 
by h cc {r), the cm. correlation function of the polymer 
system. The corresponding direct correlation function 
will be denoted c cc (r). 

The cm. correlation function h cc (r) is easily acces- 
sible in MC simulations for any polymer concentration 
IE IE 0] > but cannot be calculated directly from the 
original PRISM integral equation. However an approx- 
imate relation between h mm and h cc was established in 
Ref. which reads 

h cc (q) ~ (S^f-) h mm (q), (11) 

\»ram(?) / 

where 0, mm (q) is the previously discussed intramolecular 
monomer- monomer form factor, and VL cm (q) is the cm.- 
monomer form factor [normalized such that VL cm {q — 
0) = 1]. This relation was shown to be quite accurate 
by comparing the resulting /i cc (q), or the corresponding 
cm. structure factor S cc (q) = 1 + ph cc (q), to MC simu- 
lation data. To complete the results of Ref. [l3|. where 
only good solvent conditions were explored, and to illus- 
trate the quality of the predictions of Eq. Ill 1H for S cc (q), 
simulation results obtained for a lattice model slightly 
below its 9 temperature are reported in Fig. [S] 

However, to use Eq. ©, we explicitly need the pair 
correlation function h cc {r), i.e. we have to formulate the 
predictions of Eq. (|11|) in real space, and not in reciprocal 
space. In Fig-El we plot the simulation results of Ref. 01 
accordingly. The agreement between the "exact" h cc (r) 
and its approximation computed from Eq. i|ll|) , using the 
"exact" monomer correlation function and form factors, 
is seen to be reasonably good, except at short distances 
where the results of the approximation lie well below the 
"exact" simulation data, and can even become smaller 
than —1 at low densities, which is unphysical. Conse- 
quences of this underestimation of h cc {r) at full overlap, 
which is found under the conditions of Fig. [S] as well and 
appears to be systematic, are discussed below. 

As mentioned in Sec. [HI since there are no known an- 
alytic expressions for the two form factors appearing in 
Eq. i|ll|) for interacting polymers, we are led to use the 
expressions valid for Gaussian polymer coils in the scaling 
limit . i.e. Eq. @ and 

& cm (<l) = ^V* 2 / 12 erf (!),*= qR g , (12) 
where erf denotes the error function. In Ref. |l3j . it was 
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FIG. 5: Cm. structure factor S cc (q) for L = 500 SAW 
polymers with monomer nearest-neighbor interaction e on a 
simple cubic lattice, at (3e = 0.3 and pB? g ~ 0.04. The "exact" 
structure factor from MC simulations is compared to the pre- 
diction of the approximation (II IP using the "exact" monomer 
correlation function and form factors. 
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FIG. 6: Cm. pair correlation functions h cc (r) for L = 500 
SAW polymers on a simple cubic lattice. From bottom to 
top, pRg ~ 0.07, 0.27, 1.10. The "exact" pair correlation 
functions are compared to the predictions of the approxima- 
tion il l U using the "exact" monomer correlation functions 
and form factors. 



shown that they yield qualitatively correct results for the 
calculation of h cc (q), even in the good solvent regime. 

The importance of using exact results is best seen when 
the high density limit of h cc (q) for T > is considered, 
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FIG. 7: Zero density effective pair potentials between centers of mass (left-hand frame), for various values of the temperature 
ratio T/0; the right-hand frame shows the potentials multiplied by 4ir(r / R g ) 2 . 



which follows from Eq. 

lim ph cc (q) 



^cra (q) 
^mm (q) 



(13) 



Indeed, while this ratio is always larger than —1 if the 
exact form factors and l|12H are combined, it shows 
a region at small q where it is smaller than —1 if the 
approximation JSJ) is used. In the latter case, spurious 
singularities are thus encountered when Eq. (|SJ) is applied 
to compute c cc (q) at high densities. 

In this high density regime, an expression for the po- 
tential of mean force between the c.m.'s of polymers in 
the melt has been derived by Guenza, on the basis of 
the same thread PRISM theory as the one used here, in 
the case of excluded volume interactions only j^. Both 
approaches share similarities, like the fact that, to first 
order, the potential of mean force vanishes proportionally 
to 1/p. However, since many aspects of these theories dif- 
fer and since we are interested in the effective potential 
and not the potential of mean force, a detailed compari- 
son was not attempted. 

The outline of the present computation is thus the fol- 
lowing. First, the thread PRISM theory developed in 
Sec. [H] is used to obtain h mm (r). Then h cc (r) is derived 
by applying Eq. I|ll|) with form factors Q and (|12|) . and 
is eventually inverted using Eqs. JSJ) and © (with indices 
c instead of a) to obtain f3v cc (r). 

We first consider the zero density limit. When p = 0, 
the effective pair potential is simply the potential of mean 
force, i.e. (3v cc {r) = — ln[l + h cc (r)]. Combining Eqs. © 
(with p = 0) and (|11[1 . one finds 



f3v cc (r) = - ln[l + cRlVL cm * fl cm ( 



(14) 



where * denotes a space convolution. Results for (3v cc (r) 
and 47r(r / Rg) 2 f3v cc (r) are shown in Fig. [7\ for several 



temperatures. Above e, (3v cc (r) has a roughly Gaus- 
sian shape, with an amplitude which decreases with tem- 
perature. The potential vanishes identically for T = G 
and becomes purely attractive below, retaining a roughly 
Gaussian shape with a negative amplitude increasing in 
absolute value when T is decreased. These trends are 
globally in agreement with the behavior found in MC 
simulations jj, || , but quantitative differences are visi- 
ble. 

In the good solvent regime, for temperatures above 
approximately 11G, it is impossible to compute j3v cc {r). 
This is due to the theory producing unphysical cm. pair 
distribution functions, which are negative in the vicin- 
ity of r = 0. Moreover, approaching this existence limit 
from below, the value of (3v cc (r) at full overlap is found 
to diverge, corresponding to the fact that h cc (r — 0) goes 
to — 1 . This behavior is at odds with all available simula- 
tion data [E 0, IE IE 00 an d other theoretical approaches 
[E , showing that in good solvent the amplitude of the 
potential of mean force is always finite and has a maxi- 
mum of about 2ksT, a value reached for T ~ 50 in the 
present theory. This deficiency, which manifests itself 
throughout the high-T, low-p domain, is an obvious con- 
sequence of the underestimation of h cc (r) at full overlap, 
as illustrated in Fig. While the exact h cc (r = 0) sat- 
urates at a value strictly greater than —1 in the limit of 
infinite temperature, it is predicted by the present theory 
to go to —1.1 in the same limit. At this point it is worth 
insisting on the fact that this problem appears to be as- 
sociated with the use of Eq. pi|l. and is totally unrelated 
to the use of Gaussian chain form factors to model good 
solvent conditions. 

Around and below 0, a more complex r dependence 
of (3v cc (r) was found in the simulations of Refs. 0, ||, 
the potential globally growing more and more negative, 
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as found here, but retaining a repulsive component at 
short distances, while it is found strictly monotonous in 
the present calculation. These two types of variations 
were also found in various theoretical approaches based 
on mean-field ideas 0, |23> l34l l35l |3(j, and in fact it 
is not clear what the generic behavior of (3v cc {r) is in 
poor solvents. The available simulation data were indeed 
obtained for rather short chains, for which, on the scale 
of R g , the monomer hard cores are not of negligible size. 
This could result in the observed short range repulsive 
contribution to 0v cc (r), which would then be a finite size 
effect disappearing for longer chains |37|. 

By comparison with mean-field theories, difficulties 
can in fact be anticipated with the present theory in the 
low temperature domain. Indeed, in the zero density 
limit, following Eq. I|14fl. (3v cc (r) is a simple function of 
the overlap between the monomer distributions around 
the centers of mass of two coils at distance r. However, 
mean-field theories built exclusively on this two polymer 
overlap are known to take into account only two-body 
effective monomer- monomer interactions. This is for in- 
stance the case of the theory of Flory and Krigbaum 
Q, 13 , which can be recovered in a thread PRISM ver- 
sion by linearizing Eq. (|14fl . Such approaches are valid in 
good solvent, since then pair interactions actually domi- 
nate, but they fail in 9 and poor solvents, where inclusion 
of higher order interactions is required. By analogy, this 
raises questions on the applicability of the present theory 
in the latter domain as well. 

We now turn to the density dependence of the effec- 
tive potentials derived within the present theory. Effec- 
tive potentials obtained for various densities at T = 20 
are plotted in Fig. Q3J These results are typical of good 
solvent conditions. Starting from pR? g = 0, where the 
potential is purely repulsive and nearly Gaussian, v cc (r) 
first increases in amplitude and in range with pR g , then, 
at a density which increases when the temperature ap- 
proaches 9 from above, the trend changes: the repulsion 
at full overlap decreases as pRg is further increased, while 
the range of the effective potential still increases. In ad- 
dition, with increasing pRg, a damped oscillatory tail, 
which is more clearly seen on plots of 4Tr(r/R g ) 2 (3v cc (r), 
progressively develops when the distance between the 
c.m.'s of two coils greatly exceeds R g . All these trends are 
consistent with the simulation data of Refs. [H, IE El- 
The only qualitative difference is in the tail of the poten- 
tial, which, within the resolution of the simulations, ap- 
pears to be purely attractive. The results of the present 
theory suggest that to unveil a possible oscillatory behav- 
ior, simulations of much larger bulk systems with very 
high statistical accuracy would be required to character- 
ize the potential beyond its first minimum. 

Despite the various concerns discussed above about 
the applicability of the present theory around and be- 
low the 9 temperature (lack of more-than- two-body effec- 
tive monomer-monomer interactions, incomplete charac- 
terization of the domain of existence of the low density 
phase 30]), we present the corresponding results for com- 



pleteness. When T = 9, the effective cm. potential is 
found to vanish identically at all densities. Typical re- 
sults for the poor solvent domain, obtained for various 
densities at T = 0.759, are plotted in Fig. Starting 
from pRg = 0, where the potential is negative and nearly 
Gaussian, v cc (r) is first found to become more negative 
with increasing pRg, while roughly retaining its shape. 
Then, upon approaching the "spinodal" instability line, 
at a density which increases when the temperature ap- 
proaches 9 from below, the value of the potential at full 
overlap starts to grow with pR g and a repulsive tail pro- 
gressively develops. These tendencies are not corrobo- 
rated by the presently available simulation data. 

V. DISCUSSION AND CONCLUSION 

We have developed an integral equation approach to 
the computation of effective interactions between poly- 
mers in solution. It is based on the thread limit of the 
PRISM theory, which is used directly for the calcula- 
tion of effective interactions between average monomers, 
and in association with an approximate expression for 
the cm. pair distribution function, derived within the 
PRISM formalism as well, for the calculation of cm. -cm. 
effective potentials. 

In the case of the effective interactions between aver- 
age monomers, the predictions of the theory for poly- 
mers with excluded volume interactions only are in good 
qualitative agreement with the available simulation data. 
Moreover, we have shown that the theory lends itself to 
fully analytic calculations and compares rather favorably 
with other analytic theories like those developed for star 
polymers. 

The effective cm. -cm. interaction potentials derived 
from the present theory exhibit many of the trends ob- 
served in MC simulations. The agreement with the latter 
is very satisfactory at intermediate temperature under 
good solvent conditions, i.e. T > 9. The theory indeed 
reproduces non-trivial features like the non-monotonic 
variation of the value of the potential at full overlap 
and the development of a long range tail when the den- 
sity increases. This agreement deteriorates at very high 
temperature, in particular for (3 = (SAW model), and 
around and below the ^-temperature. 

Since the HNC closure used to compute the effective 
pair potential is exact for all practical purposes, when 
applied to "soft-core" systems which exhibit essentially 
"mean- field fluid" behavior '27], the observed discrep- 
ancies must be traced back to the other approximations 
made in the theory, namely the RMPY/HTA closure used 
to compute h mm (r) in the thread limit, the use of form 
factors Qmm{q) and Q cm (q) valid for ideal rather than in- 
teracting chains, and Eq. Ijllll relating h cc (q) to h mm (q). 
Simulation data tend to show that the use of Eq. I jllj l is 
the main cause of the breakdown of the theory at very 
high temperature, while there are indications that the 
present use of the RMPY/HTA closure in the low tern- 
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FIG. 8: Effective pair potentials between centers of mass (left-hand frame), at temperature T = 20, for various values of the 
reduced density pR g ; the right-hand frame shows the potentials multiplied by An(r / R g ) 2 , and the insert shows an enlargement 
of this function for large r/R g . 
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FIG. 9: Effective pair potentials between centers of mass (left-hand frame), at temperature T = 0.750, for various values of 
the reduced density pR 3 g \ the right-hand frame shows the potentials multiplied by 4ir(r/R g ) 2 . 



perature regime is not justified. In the present state of 
the theory, the use of ideal chain form factors is in fact 
not found to be a critical approximation. 

In conclusion, we have proposed a first step towards 
an integral equation theory of effective pair potentials 
between macromolecules in solution, which could bypass 
costly simulations for the determination of the cm. pair 
correlation function. The results appear very encourag- 
ing, since, while only the simplest ingredients have been 
combined in the present theory, and there is clearly room 
for improvement, it is already found to be able to cap- 
ture various non-trivial behaviors observed in recent sim- 



ulations. It appears thus that integral equation theories 
could be in the future a valuable tool for the development 
of coarse-grained descriptions of complex polymeric sys- 
tems. 
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APPENDIX A: THREAD PRISM CALCULATION 

In this Appendix, we present our calculation of the 
monomer pair correlation function based on PRISM the- 
ory in the thread limit. It is inspired by the work of 
Chatterjee and Schweizer 0>E3 (see also Ref. ^j) for a 
related approach), a special asymptotic regime of which 
is obtained. We start with Eq. J3J, involving quantities 
which all remain finite in the thread limit as defined in 
Sec. [n] The case of excluded volume interactions only is 
considered first, and then used as a reference system to 
include attractive interactions between the monomers in 
the closure relations. 



1. Reference system: Excluded volume interactions 

Considering first the excluded volume interaction only, 
a Percus-Yevick (PY)-like closure is used, in which 
C mm (r) — c R g S(r) (R g is introduced here to make cq 



dimensionless) , or equivalently in Fourier space, 



(Al) 



This singular C mm (r) is required to account for the effect 
of excluded volume on intermolecular correlations in the 
limit of vanishing size of the hard core of the monomers 
in the thread limit. 

The PRISM-OZ relation © then leads to Eq. ®, with 
c = Co, and the value of Co is determined by enforcing the 
hard core condition h mm (r = 0) = — 1. 

In general, all the calculations have to be performed 
numerically. This is the case in the present work when 
!)„((]) is chosen to be the Debye function Q. However, 
fully analytic results can be obtained if the usual Pade 
approximant of this function, Eq. J5J, is used. Then, one 
finds 



c = -27iV2 - 2tt 2 P R 3 



and 



e -r/t„ _ e -r/i c 



r/R g 



(A2) 



(A3) 



2irpR 3 g 

where £ c = R g /V2 and £ p = R g /(V2 + 2wpR 3 g ). 

2. Including an attractive tail 

We now introduce the attractive interaction between 
the monomers, which, for convenience, is taken of the 
Yukawa form QJ. In Ref. 15], various closures of the 
PRISM-OZ equation have been studied in the thread 
limit, for chains with monomers interacting through a 
hard core and an attractive tail. Here, we consider the 
reference molecular Percus-Yevick /high temperature ap- 
proximation (RMPY/HTA) El, wnicn has been 



found to be the most satisfactory in this earlier work. 
With the present notations, it reads 



C mm (r) = C re f(r) - (3V{r)g rct (r), 



(A4) 



where a renormalized attractive potential V(r) = L 2 v{r) 
was introduced. C re f(r) and g-cefij') are respectively the 
renormalized monomer direct correlation function and 
the monomer pair distribution function of the reference 
system, the calculation of which is explained above. 

As mentioned in Sec. [HI on physical grounds, we sup- 
plement the previously applied limiting process with the 
limit of a vanishing range of the attractive potential, 
namely a — > 0, keeping the ratio A = a/a constant. 
To take this supplementary limit in the RMPY/HTA 
closure, we need to evaluate \im a ^, (3V(r)g rc f(r). We 
assume that g r ef(?") satisfies the hard-core condition, 
<?rcf(0) = 0, and has a finite slope at the origin. This 
hypothesis holds both for the analytic result, Eq. i|A3() . 
and when the exact Debye form factor is used instead of 
its Pade approximant (see below). Then we can write 
the Taylor expansion 

g ie{ (r) = g' iei (0)r + o(r). (A5) 

Therefore 

f3V(r)g ie{ (r) = -L 2 0eae- r / a g' iei (O) + o(l) 

e -r/a (A6) 
= -36/?eA 4 i? 4 g > (o) + (i), 

where R 2 — La 2 / 6 has been used. Now, taking the limit 
a — > 0, we find 

lim (3V(r)g Tei (r) = -288^eA 4 i?V ef (0)<5(r), (A7) 

since the contribution of the o(l) in (|A6I) goes to zero. 
The RMPY/HTA closure relation hence reduces to 



C m m( r ) — c o 



l-|)/OMr). 



with the temperature O given by 



9 = -288ttA 



4^(0)^ 



(A8) 



(A9) 



[remember that Co < 0, as seen in Eq. i|A2|l ]. The latter 
expression can be greatly simplified, using the fact that, 
for large q, h lc f(q) decays like 4co/ '(g 4 i? g ); this allows to 
show that 



<4f(o) = - 



C() 



2nR 



which yields 



6 = 144A" 



(A10) 



(All) 



Thus, we find that is the same for both choices of the 
polymer form factors, Eq. (0J or (jSJ), and is independent 
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of the polymer density. Fourier transforming C mm (r), 
defining c as in Eq. Q, and injecting the result in Eq. 
one obtains Eq. ®. 

As anticipated by its notation, is in fact the 9 tem- 
perature of the present polymer model. This can be 
first seen from Eq. I|A8|I . since one finds that C mm (r), 
often interpreted as a renormalized two-body monomer- 
monomer interaction, vanishes for T = 0. More evi- 
dently, this can be shown by considering the pressure of 
the polymer fluid. Following the "compressibility route" 
pfij . the pressure is indeed easily calculated to be 



performed, leads to 



P- 



nV2R 3 g 



^ 2 I?6 



(A14) 

which is precisely the result obtained by combining 
Eqs. ljA12|) and (|A2|) . This indicates that the present 
theory corresponds to an a — > limit of the theory of 
Chatter jee and Schweizer. 



Limitations of the model 



PP= / [1- p'C mm (q = 0,p')]dp' 

rP (A12) 



p'c (p')dp', 



and one sees that, for T = 0, the pressure is that of a 
perfect gas, meaning in particular that its second virial 
coefficient vanishes, which is the usual definition of the 6 
temperature. 

Here again, a fully analytic expression for h mm {r) can 
be obtained if the Pade approximant of the chain form 
factor is used. It retains the form i|A3(l . with £ c un- 
changed and £ yP = R g I J 2(1 — pcR?). 



Relation to the theory of Chatterjee and 
Schweizer 



We show here that the present calculation is related 
to an asymptotic regime of the theory developed in 
Refs. [lj, llfij . in which the range of the attractive tail 
was kept finite. To this purpose, the calculation of the 
pressure is considered. 

In Refs. |l5l Il6j |. where the approximate chain form 
factor (0 was used, the pressure was found to be given 
by (see Eq. (30) of Ref. p3) 



PP = p + irV2R 3 Q 



J 2 E>6 „3 

2 R aP 



flex 



( 127rV6A 3 i^VZp 2 W&LP 

\ 9 l+(aV2/l 



VqalVl 



nRl 



In 



V2/R g ) 
2nV6AR 3 g p 
L[l + (ay/2/R g )] 



(A13) 



where the notations of the present paper are used and 
a and p m have been reexpressed as functions of L, R g , 
and p, using the relations R 2 g — La 2 /6 and p m = pL. 
a can be handled in a similar way, using the fact that 
R 2 = La 2 /(6A 2 ), and all the variables which vanish in 
the thread limit as defined in Sec.|n]are then eliminated. 
The only remaining singular quantity is L, which should 
go to infinity. Taking this limit in Eq. I|A13(1 . after a 
large L expansion of the logarithm and denominators is 



The present polymer theory shows some limitations 
which are reviewed here. 

First, in the presence of attractive interactions, the 
hard core condition h mm (r = 0) = — 1 is no longer en- 
forced 15]. For instance, one finds for the analytic solu- 
tion 



Ar = 0) = 



1-Jl-pcR* 



irV2pR g 



(A15) 



which is strictly larger than —1 for all densities and fi- 
nite temperatures. Such a feature is very common in this 
kind of thermodynamic perturbation type approxima- 
tions, where the direct correlation function is completely 
specified for all values of r. We expect its impact to be 
modest on the outcome of the proposed coarse-graining 
procedure, since we will be interested in a description on 
a mesoscopic level, while this inaccuracy manifests itself 
on the microscopic, monomer level. 

Secondly, below the 6 temperature, there is a domain 
where the integral equations exhibit a singularity, corre- 
sponding to the vanishing of the denominator of Eq. I©. 
This is in fact the remnant, when L goes to infinity while 
taking the thread limit used in this work, of the finite L 
spinodal instability found by Chatterjee and Schweizer 
[l5j, with a critical temperature equal to and an in- 
finite critical value of pR g . The boundary of this low- 
temperature, high-density domain can be explicitly com- 
puted when the approximate expression of the form fac- 
tor is used, leading to 



pRl 



1 



ttV2 







e-T 



- 1 



T < 0. 



(A16) 



This means that, below the 9 temperature, where poly- 
mer solutions are known to phase separate between 
polymer-rich and solvent-rich fluids, the proposed theory 
will only allow to investigate the latter phase. 



APPENDIX B: ANALYTIC EFFECTIVE PAIR 
POTENTIAL BETWEEN AVERAGE 
MONOMERS 

A fully analytic expression for (3v aa (r) can be obtained 
if the approximate expression (J5J for the chain form fac- 
tor is used. Indeed, the expression for h aa (r) = h mm (r) 
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is known; solving Eq. JSJ for h aa (q) and identifying the 
result with Eq. 10, one finds 



If c > 0, i.e. T < 8, c aa (r) shows damped oscillations, 



c aa (q) 



1 - pcRffi mm (q) + pcR 3 Sl mm (q) 2 



( B1 ) c aa (r) = - 



Replacing fl mm (q) by its expression JSJ), c aa (q), and sub- 
sequently c aa (r), are easily obtained. Two regimes are 
found, depending on the sign of c. 

If c < 0, i.e. T > 0, one finds that c aa (r) is the 
difference between two Yukawa terms, 



Caa(r) 



2nJpcm(pcm-4) 



p — »7£+ _ p —r/t- 

x - jw ' ( B2 ) 

r/Rg 



xl^lsin(r/6), (B4) 



nJpcRUA-pcRl) r / R 9 



where 



& = R 9\l 1 — ~~Bs' & = R 9il-=L- (B5) 



4 - pcRjj ' 



where 



Combining these results with Eqs. (| A3|> and © pro- 
(B3) vides the desired analytic expression for (3v aa (r). 



2-pcRlTJpcRl{ P cRl-A) 
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